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ANALYSIS OF ARCHES BY FINITE DIFFERENCES 


Ephraim G. Hirsch,’ J.M. ASCE, and E. P. Popov,? A.M. ASCE 


SYNOPSIS 


In flexible arches vertical and horizontal displacements of the rib due to 
imposed loading cause large increases in the bending moments computed on 
the basis of invariant arch geometry. These increases in moments cannot be 
neglected in safe design. Present methods for solving such problems are 
either restrictive or are cumbersome in application. The present paper is 
an attempt to remedy this situation. Thus, starting with two basic differential 
equations which relate both the vertical and horizontal deflections of the rib to 
the internal forces, a solution of these equations by the finite difference 
method is discussed. Arches of any geometry, cross-section, and conditions 
of loading may be analyzed by the proposed method. Moreover, the proposed 
method is rapid in applications, and, due to its simplicity, may be mastered 
readily. 


NOTATION 
The following list contains the most important of the symbols used in this 
paper. 
cross-sectional area 
modulus of elasticity 
maximum rise of arch 
horizontal thrust 
length of finite difference interval 
moment of inertia 
span length =nh 
bending moment of arch rib 
simple beam bending moment 
fixed end moment 
bending moment due to horizontal displacement 
bending moment due to vertical displacement = Hn 
axial thrust 
number of finite difference intervals 


1. Structural Designer, I. Thompson, Cons. Engr., San Francisco, Calif. 
2. Prof. of Civ. Eng., Univ. of California, Berkeley, Calif. 
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t change of temperature 

Vv R left, right vertical reactions 

w(x) transverse loading 


dead load, load per unit length 
wf live load, load per unit length 
x,y Cartesian co-ordinates defining original shape of arch axis 
a coefficient of thermal expansion 

vertical deflection of arch axis 
slope of arch axis 
horizontal deflection of arch axis 


INTRODUCTION 


There are two basic theories that may be employed in the analysis of arch 
bridges. The classical theory of analysis, often known as the “elastic” theory, 
assumes that any deflections that occur in the rib due to external loading are 
so small as to have negligible effect upon the internal bending moments. 

Such an approach is suited only to the analysis of very rigid masonry or con- 
crete arches. However, the design of flexible arches requires consideration 
of the changes in geometry of the rib that occur when external loadings are 
imposed upon the structure. It is no longer possible in these cases to con- 
sider that the arch has an invariant geometry, for the changes in geometry, 
though small in comparison with the co-ordinates of the arch, in some cases 
cause bending moments that are much larger than those computed by elastic 
theory. 

Consider the two-hinged arch shown in Fig. 1 which is acted upon by some 
arbitrary vertical loading w(x). The resultant reactions at the supports will 
be Vp and Vi, vertically and a horizontal thrust H. A point A on the rib will 
displace 7 in a vertical direction and & in a horizontal direction to a new 
position A'. The co-ordinates of point A are x and y in the undistorted posi- 
tion of the arch and x + — and y + 7 in the distorted position A’. In the elas- 
tic theory it is considered that the bending moment at any point on the arch 
rib is the difference between the bending moment that would exist if the span 
were a simple beam and the moment caused by the product of the horizontal 
thrust and the initial ordinate of the point in question. However, it can be seen 
that any horizontal deflections will change the lever arms of the vertical loads 
and thus will change the value of the simple beam bending moment. More- 
over, any vertical deflection will change the lever arm of the horizontal thrust. 
Now, while the numerical values of the simple beam moment and the moment 
caused by the horizontal thrust are in themselves quite large, the magnitude 
of their difference becomes of the same order as the products Hy, Vé , and 
wix)é . Small deflections, therefore, can drastically change the net design 
moments. In such cases it is necessary to employ an analysis that will in- 
clude the effects of deflection in the solution of the problem. This type of 
analysis is known as the “deflection” theory of arch bridges. For the correct 
and safe design of large span, flexible arches, the deflection theory is indis- 
pensible and it does not represent a mere refinement of the elastic theory. 
Initial investigations of the deflection theory of arches were made ina 
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series of papers(!,2,3,4,5)8 published by Josef Melan and Fritz Engesser be- 
tween 1900 and 1906. In 1934 B. Fritz‘) extended Melan’s work to cover 
three-hinged, two-hinged or hingeless arches. In a theory developed in 1931, 
s. Kasarnowsky\7) presented a simple calculational procedure. In 1935, A. 
Freudenthal (8,9) proposed a theory for a two-hinged or hingeless arch in 
which the variable moment of inertia I was equal to I, sec @ , where I¢ is 
moment of inertia at the crown. 

In all of the above procedures the horizontal component of deflection is 
neglected. This is a satisfactory assumption only for very flat arches. In 
addition, all of the enumerated procedures are restricted to arches which 
have a definitely prescribed variation of the cross-section or, except for the 
work of Fritz, have restrictions as to the placement of live loading upon the 
structure. To remedy the situation, G. Stern(10) developed a method utilizing 
a Fourier series solution for a two-hinged arch and was able to solve the 
problem for any variation of arch rib or loading. However, the horizontal 
component of deflection was still neglected. This work was extended by W. F. 
Lu(11) in 1951 to a hingeless arch. A simple approximate procedure by use 
of a magnification factor was developed by Hardesty et al(12) in conjunction 
with the Rainbow Arch design in 1945. The same point of view was also 
used(13) by R. S. Rowe. 

In 1948, P. N. Chatterjee(14) developed a method for the analysis of two- 
hinged arches by use of the Newmark numerical procedure. The advantages 
of his work are that it may be applied to an arch with any variation of rib 
properties and live loading and, in addition, it accounts for the horizontal 
component of deflection. This work was checked by means of a model study 
by C. L. Monismith(15) in 1954. 

Any attempt to solve the deflection theory of arches for a general case of 
arch rib properties and live loading presents an almost intractable mathemati- 
cal problem unless some sort of numerical procedure is employed. One such 
approach is Chatterjee’s solution, cited above. It is, however, a successive 
approximation type solution and requires several cycles of computation to ob- 
tain the final results. It is proposed in this paper to demonstrate that the use 
of the finite difference method of solving differential equations allows the rapid 
and accurate analysis of any arch structure regardless of its geometry, cross- 
sectional variation or external loading. It is furthermore applicable to either 
the elastic theory or deflection theory of arch action. The special steps and 
precautions required in the application of this method to the arch will be dis- 
cussed first, and then, in order to aid in clarification of this discussion, 
numerical problems for both the elastic and deflection theories will be pre- 
sented. 


Differential Equations of the Arch 


The action of arch ribs is governed by two differential equations which re- 
late the vertical and horizontal deflections of the rib to the internal forces in 
the rib. The two equations are not widely known and to the best of the authors’ 
knowledge, have not been cited previously in civil engineering literature. 
These equations were developed(16) and are stated here in terms of the hori- 
zontal and vertical deflections, rather than in terms of radial displacements. 
The two equations are: 


3. Numbers in parentheses refer to the Bibliography at the end of the paper. 
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The first term on the right-hand side of Eq. 1 is due to flexure while the 
second term is due to axial deformation. If flexural effects alone are to be 
considered, then this latter term is omitted. Note that if the effects of axial 
deformation are neglected and thus only the first term on the right-hand side 
of Eq. 2 is retained, the equation becomes one which relates the rate of change 
of vertical deflection to rate of change of horizontal deflection. This equation 
is then in effect a “coupling” equation. 


Solution of the Arch by Finite Differences 


. The Basic Procedure 


For any given arch and condition of loading, it is necessary to solve the 
basic Eqs. 1 and 2 for the deflection functions & (x) and n(x) in order to deter- 
mine the distorted geometry of the arch. In addition, the redundant reactions 
must be determined. It is then possible, as may be seen from Fig. 1, to 
evaluate by statics the internal forces at any section of the arch rib. 

In this paper it is suggested that the solution be accomplished by the use of 
finite differences. The characteristic of this procedure is that the differen- 
tial equations are transformed into a series of algebraic equations. To ac- 
complish this, approximate relations for the derivatives are used. For exam- 
ple, in the “three-point-central difference” approximation, the first and 
second derivatives of a function f(x) at an m-th point may be expressed in 
terms of the values of the function at the adjoining points as follows: (17 


To express higher order derivatives it is necessary to use five or more 
adjoining points because of the restrictions inherent in the derivation of Eqs. 
3 and 4. However, since the arch differential equations are not higher than 
second order, the three-point approximation is sufficient. 

To apply the method of finite differences, the span of the structure to be 
analyzed is divided into any desired numbe~ of segments. For ease of solu- 
tion, it is usually advisable to have segments of equal length, although a mesh 
of variable spacing may be used to advantage if the behavior of a particular 
point on the arch needs to be investigated more accurately. In any event, the 
degree of accuracy of the solution is dependent upon the fineness of the mesh 
chosen, Excellent results were obtained for the numerical examples of this 
paper by using a mesh of only 8 subdivisions. 

At each of the interior points obtained by this subdivision the differential 
equations are converted to algebraic equations by using the finite difference 
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expressions for first and second derivatives and by substituting the values at 
each point of the physical constants occurring in these equations. On this 
basis, if the arch span is divided into n segments, (n - 1) algebraic equa- 
tions are obtained for each of the basic differential equations. However, 

the set of algebraic equations formed from Eq. 1 will contain n unknowns: 

an unknown vertical deflection, 7, at each interior point, and the hori- 
zontal thrust, H. It is, therefore, necessary to eliminate one unknown from 
this set of algebraic equations by the use of Eq. 2 and the boundary condi- 


tions of the problem. 
The Boundary Conditions 


For a pin-ended arch, since the moments at the ends are zero, the bound- 
ary conditions in finite difference form become: 


1) 3 at left springing 


(5) 


where the subscript -1 signifies an imaginary point one interval h to the left 
of the left hinge. The subscript n+1 signifies an analogous point to the right 


of the right end. 
For a hingeless arch, the boundary conditions are: 


(6) 


n+ 


These boundary conditions may be used in the case of a hingeless arch to 
define the redundant end moments. Thus, writing Eq. 1 in finite difference 
form at the left springing and neglecting axial distortions, the following rela- 
tion is obtained: 


aul %,=O, and 


sec 

A similar expression may be formulated for the moment at the right spring- 
ing of an arch and these expressions may be used in lieu of the redundant un- 
known moments in the solution of a problem. 

When Eq. 2 is applied in finite difference form at the boundary points in 
addition to the interior points of the subdivision, it becomes possible to ex- 
press vertical displacement in terms of each other. This permits the elimina- 
tion of one of the superfluous n's from the algebraic equations based on Eq. 1. 
In order to avoid a trivial relationship of zero equals zero, for hingeless 


arches it is first necessary to differentiate Eq. 2 once with respect to x before 
writing it in finite difference form at the boundary points, In addition, Eq. 2 
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may be used to define all horizontal displacements in terms of the vertical 
displacements. This is necessary to eliminate the terms that appear in 
Eq. 1 when horizontal displacements are considered in the analysis. 


Procedure for the Elastic Theory 


In elastic theory the object of the analysis is the determination of the re- 
dundant forces. The analysis is not concerned with the deflections of the arch 
per se. Although it recognizes that the deflections exist, it assumes that they 
have no effect upon the final stresses. For a two-hinged arch, the algebraic 
equations, the formulation of which was discussed above, can be solved for the 
value of the horizontal thrust H. For a hingeless arch, it is also necessary 
to determine the vertical deflections occurring at each of the stations adjacent 
to the two springings. These deflections may then be used to determine the 
two redundant end moments by means of Eq. 7. 


Procedure for the Deflection Theory 


If an arch is to be analysed by the deflection theory, the expressions for 
internal bending moment must first be modified to account for the change in 
geometry of the rib, It was indicated in the Introduction that vertical deflec- 
tions affect the moment caused by horizontal thrust and horizontal deflections 
change the moments caused by vertical loading. The lever arm for the thrust 
is modified at each station by expressing it as the ordinate, y, at the station 
plus the as yet unknown vertical deflection, 7 , at that station. In algebraic 
form, this changes the expression for bending moment to: 


M-M'-H(¢ +7) 


No such general rule is possible when dealing with the horizontal displace- 
ments. For general cases of transverse loading, it is necessary to formulate 
the primary bending moment at each point by statics, considering that each 
point has a positive horizontal displacement —. For the special case of an 
arch carrying a uniformly distributed vertical load, however, it is possible to 
obtain a reasonably accurate general expression for the primary bending mo- 
ment modified for the change in geometry. Thus, consider the arch shown in 
Fig. 1 which carries a partial uniform vertical load w causing the vertical re- 
actions V;, and Vp. These reactions will be assumed to be those for a simple 
beam carrying the indicated transverse load. The simple beam moment about 
A‘ will then be 


M's Vi Z +B) 


§,(V,- 


This expression holds for situations where the uniform load extends from the 
springing at least to the point in question. A similar expression can be formu- 
lated if the load extends from the right-hand springing. 
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Note that the consideration of vertical deflections gives rise to terms 
which are the products of the horizontal thrust and the vertical deflection at 
each point. Both these quantities are unknown and it is desirable that one of 
them be intelligently estimated in order to permit the direct solution of the 
matrix of the algebraic equations based on Eq. 1. The effect that deflections 
have upon the value of horizontal thrust is usually relatively small. Hence, 
the arch may be solved initially for the value of this redundant force by the 
elastic theory and the value of this thrust then be used in the deflection theory 
analysis. 

The above modifications in the expression for the bending moments are all 
that are required for the deflection theory solution. If horizontal deflections 
are not to be considered, then the change in the lever arm of the horizontal 
thrust alone need be considered. In addition, since the thrust is assigned an 
assumed numerical value, the only unknowns that then appear are the various 
vertical deflections. If horizontal deflections are considered, it is necessary 
to eliminate the additional unknowns by the methods described previously. 

Upon solving the algebraic equations for the various deflections, the bend- 
ing moments in the arch may be determined by the substitution of these values 
into the previously formulated Eqs. 8 and 9. The true bending moments are 
thus found directly from the initial solution and several successive approxima- 
tions are not required. However, initially the value of thrust based on the 
elastic theory was assumed and does not exactly correspond with its true 
value. It is likely that any difference between the two values will be insignifi- 
cant. If it is desired to solve the equations by an iterative procedure rather 
than by the direct solution suggested here, the H may be left initially as an un- 
known. Such a procedure may be carried out with the aid of electronic 
computers. 

It should be noted that if for symmetrical loading, the *three-point central 
difference” approximation is used in establishing the finite difference equa- 
tions the apparent result is that equal horizontal displacements occur for the 
first point on each side of the crown. This contradicts the actual situation of 
equal and opposite horizontal displacements at these two points. This con- 
tradiction can be eliminated by using the “five-point central difference” ap- 
proximation in setting up the finite difference equations. Alternatively, with 
the “three-point” approximation the finite difference equations may be written 
from both ends of the arch toward the crown. Then a relation for the horizon- 
tal displacements at the crown point need not be used at all and no contradic- 
tion occurs, 


Numerical Examples 


The two-hinged arch shown in Fig. 2 will be solved in this section by both 
the elastic and deflection theories, In both solutions the effects of axial dis- 
tortion will be neglected in order to obtain values which can be compared with 
other known analyses. For this purpose, the differential Eqs. 1 and 2 simplify 
to 


(1a) 


MM! 
and 

LF =_— Aq C (2a) 
AX AX 
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Elastic Theory 


To solve this problem by the finite difference method, the arch is divided 
into 4 segments. The geometric properties of the arch, including the primary 
bending moment, at the various points are tabulated in the first six columns 
of Table 1. Since deflection effects are neglected in this type of solution, 

Eq. 8 becomes M = M' - Hy. 


TABLE 1 


1 


2 2 
n= kh; h = 2 ft; BI = 1.832 lb-rt*; = 2.18 x 107 1b” 


Mt 


tan @ sec @ (ft-1b) 


Sta. 


0.6250 1.1792 
0.9375 0.3125 1.0478 56.10 
1.2500  0,0000 1.0000 691.20 
0.9375  =0.3125 1.0478 88.10 


0 0.6250 1.1792 0 


Since this arch has pinned ends, the boundary conditions are given by Eqs. 
5, and Eq. 1 in finite difference form may be written as follows: 


At station 1: 


-21,+1%,) = sec (1, - Hy,) 


and since, 0, 


( “21, + 1, ) = 2./8x/0 |, 0478/. 549,40- 0.93154) 


or 


(/254.6 -2,/448H)/0” 


At Station 2: 


(%, Ns) = 2,1BX/0 x 1.000 (69.20-/ 250H) 


-2g + No, = (/50%2-2.7293 H) 
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M 
0 
1 30.0 
22-5 
2 
=-22.5 
3 =30,.0 
=22.5 
|| 
or 


At Station 3: 


(11, “27s + %4)= 2.18 %/0 %/0978 [488.40 


and since n4 = 0, 


There are four unknowns 11, 2,13, and H, but only three equations. It, 
therefore, is necessary to utilize Eq. 2 together with the boundary conditions 
to obtain a solution for H. Writing Eq. 2 in finite difference form, and then 
using Eqs. 5, the following relations are obtained: 


At Station 0: 


= - 9. 6257, 


At Station 2: 


At Station 4: 


(6-€,)=-( -0.625) ov = 0.625, 
and since > Y%4=-F,- 


This condition allows the elimination of one unknown from the group of rela- 
tions, Eqs. 10, formulated from Eq. 1. There are now as many equations as 
there are unknowns and they may be readily solved for the result H = 552.96 lb. 
It can be shown that this value is exact, even though a very coarse mesh was 
used, 


Deflection Theory 


In this solution the effects of both vertical and horizontal deflections will 
be included. Accordingly, Eqs. 8 and 9 must be used to formulate the bending 
moment at every point; the remainder of the solution is identical with the 
elastic theory soiution. However, now let the arch be divided into eight equal 
segments of length h = 1 ft. The geometric properties of the arch, together 
with the primary bending moments modified for horizontal deflection effects, 
are given in Table 2. 

Before applying Eq. 1, it is advisable to utilize Eq. 2 to obtain expressions 
for in terms of 7. Thus, utilizing Eqs. 5, the following relations are ob- 
tained when Eq. 2 is converted to finite difference form: 


At Station 0: a(é- -E )=-3 + (4 - 1.,) ton 8, 
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Sta. 


0 
1 
2 


and since 


At Station 2: 


and since 


or 


6;he1 ft; 


At Station 1: 


At Station 3: 


x y 
(ft) (ft) 


0.Sh69 
0.9375 
1.1719 
1.2500 


1.1719 
0.9375 
0.569 


an wn & 


7. 


( 


EI = 1.632 1b-rt*; 


2 


tan 9 sec 


0.62500 1.1792 
1.1015 
0.31250 1.0478 
0.15625 1.0122 
0 1.000 
0.15625 1.0122 
-0.31250 1.0478 
l:6875 1.10h5 


E =- 0.625% 


and “77, =O, 
=- 0.46875 


=-0.6257, , 
=- 0, 3125 (%45+7%,) 


(%4-2)* 0./S625 (12) 


(14424) 


at 
fr” 0.5585 x 10 ~ lb 


-1 


0 
324290 274620 5, 
5L8.L0 + 172.606 
670.50 
691.20 = 30.00 54 
625.50 = 101.40 
488.0 172.80 
279290 = 


iit 


At Station 8: 57) =-(1,- 7) x (- 0.625) 


= 0.625%, 


At Station 7: (€, (Ne-%e)*(- 0.46875) 


and since oO and O, 


= 0.46875 
At Station 6: (E,-E5)=-(%-% x (-0, 3/25) 


0.3/25 (4, +7.) 


Eq. 1 may now be written in finite difference form at the seven interior 
stations. With the-aid of Eqs. 12 this will result in seven equations with 
seven unknowns, provided it is assumed that H = 552.96 lb. This value of H 
is based on the elastic theory solution, An example of the conversion of Eq. 1 
to an algebraic equation is made below for station 1: 


0.545 85410 1.10 45] (324.90 + 274.20 ) 
- 552.96 (0.5469+7,) | 


But since from the first of Eqs. 12 € , = -0.625n, , the above equation reduces 
to 


-—/563,303 /000 Te = /3,556 (13) 


For the remaining six interior points, the procedure for setting up the al- 
gebraic equations is similar. The final matrix of coefficients of these equa- 
tions is shown in Table 3. For this matrix Gauss’ Scheme, also known as 
Doolittle’s Method, was used to obtain the solution of the simultaneous equa- 
tions. The large number of zero terms lying outside of the main diagonal 
greatly expedite the solution. The unknown £ terms are determined from 
Eqs. 12. 

The results are assembled in Table 4 and are plotted in Figs. 3, 4, and 5. 
The values in the last column of Table 4, the internal bending moments, are 
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obtained from a sum of the primary bending moment, the bending moment 
contributed by horizontal deflection, and the bending moment caused by the 
horizontal thrust. For example, at Station 1: 


M- M'+M¢ -H(y4 +7) 
324.90 + 274,20x0.037/ 
552.96(0,5469- 0.0626) 


324,90+/0,72 -267,90 
+67.82 ft-/b 


Table 2 
Matrix of Coefticierts for Equation | 
Leflection Theory 


h 
1s G Ey Msec9 
4 
0 


[503.303 1000 


1000 1/000 0 

+/2.3Z8 1000 -l68Z/56 1/000 0 
QO 1000 -7007Z5 1000 
0 1000 -Ié7e.976 1000 +/75086 +/Z439 
0 1000 /000 +/7:/59 
O 0 /000 -/$74.607| \+/2574 


-/2,.474 


TABLE 


M 


(ft) (ft-lb) 
0.0391 6762 
0.0395 63 elu 
0.0378 57 
0.0264 - 0.51 
0.0371 57 67h 
0.0387 ~82,31 
0.0381 
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a? 
7s 
Stae 

1 0.0626 
2 0.0813 
3 0.0583 
0.0005 
5 0.0569 
6 00.0825 
7 0.0610 

ae 


Also plotted on Figs. 3, 4, and 5 cre the results of a finite difference solu- 
tion similar to the above except that mesh of 16 was used. In addition, the 
results from an analysis using the Newmark procedure (15) are also shown. 

It is significant to note that the analysis with a coarse mesh of only 8 sub- 
divisions yields values that are in such close agreement with the answers 
from more accurate solutions. At the point of maximum moment, there is 
only a 6 per cent discrepancy between the values from the 8 and the 16 divi- 
sion solutions. However, as may be seen from Fig. 5, the bending moment 
curve based on the elastic theory is considerably in error. For example, the 
bending moments at the quarter points are over 100 per cenit inerror. Quali- 
tatively, this may be an extreme case, but by no means an impossible one. 

The fact that the effects of axial strain were neglected does not seem to 
have altered the accuracy of the final answers, as may be judged by a com- 
parison of analytical and experimental results. As is known, axial strains 
effects become of importance in the case of the very flat arches. Their in- 
clusion in the finite difference solution presents no serious difficulties. An 
examination of Eqs. 1 and 2 shows that, with the exception of the term for 
axial thrust, N, merely numerical constants are contributed to the right hand 
side of the corresponding equations written in finite difference form. The 
term N is a function of the vertical shear and horizontal thrust at any point 
and thus may be expressed in terms of a numerical constant and the unknown 
thrust H at each point. 


CONCLUSIONS 


As indicated by the foregoing discussion, the application of finite differ- 
ences to the differential equations that govern the action of arch ribs makes 
possible a rapid and accurate solution of these equations. The solution is a 
completely general one which is not restricted in any manner by variations 
in external loading or geometrical properties of the arch itself. The degree 
of accuracy of the results is determined by the fineness of the mesh chosen. 
It has been shown that suprisingly good results were obtained even with the 
use of mesh of 8 sub-divisions. The procedure may be used in conjunction 
with either the elastic or deflection theory of arches. This is dependent upon 
whether or not the expression for the bending moment is modified initially for 
effects of horizontal and/or vertical deflections. If the expression is not so 
modified then solution of the differential equations by the finite difference 
method is an alternate to the usual methods of determining the redundant 
forces of indeterminate arches, Actually, in some cases it may be much 
simpler to apply this method than the customary procedures, 

However, it is felt that the finite difference method has its greatest applica- 
tion in the solution of the deflection theory problem. Here, it definitely pos- 
sesses the advantages of rapidity, accuracy and generality over the other 
available methods. With the advent of the electronic and digital computers 
for the large scale inversion of matrices, the finite difference method should 
prove particularly advantageous for solutions of large arch structures. 
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